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Abstract. We derive local integral and sup- estimates for the curvature 
of stably marginally outer trapped surfaces in a sliced space-time. The 
estimates bound the shear of a marginally outer trapped surface in terms 
of the intrinsic and extrinsic curvature of a slice containing the surface. 
These estimates are well adapted to situations of physical interest, such as 
dynamical horizons. 


1 Introduction 

The celebrated regularity result for stable minimal surfaces, due to Schoen, Si¬ 
mon, and Yau [SSY75], gives a bound on the second fundamental form in terms 
of ambient curvature and area of the surface. The proof of the main result of 
[SSY75] makes use of the Simons formula [Sim68] for the Laplacian of the sec¬ 
ond fundamental form, together with the non-negativity of the second variation 
of area. In this paper we will prove a generalization of the regularity result of 
Schoen, Simon, and Yau to the natural analogue of stable minimal surfaces in the 
context of Lorentz geometry, stable marginally trapped surfaces. In this case, a 
generalization of the Simons formula holds for the null second fundamental form, 
and the appropriate notion of stability is that of stably outermost in the sense of 
[AMS05, New87]. A local area estimate for stable marginally trapped surfaces, a 
generalization of a result due to Pogorelov [Pog81] allows us to give a curvature 
bound independent of assumptions on the area of the surface. 

*Supported in part by the NSF, under contract no. DMS 0407732 with the University of 
Miami. 
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Let E be a spacelike surface of co- dimension two in a 3+1 dimensional Lorentz 
manifold L and let be the two independent future directed null sections of the 
normal bundle of E, with corresponding mean curvatures, or null expansions, 9 ± . 
E is called trapped if the future directed null rays starting at E converge, i.e. 
9 ± < 0. If L contains a trapped surface and satisfies certain causal conditions, 
then, if in addition, the null energy condition is satisfied, L is future causally 
incomplete [Pen65]. Let / + be the outgoing null normal. If L is an asymptotically 
flat spacetime this notion is well defined, otherwise the outgoing direction can be 
fixed by convention. We call E a marginally outer trapped surface (MOTS) if 
the outgoing lightrays are marginally converging, i.e. if 6 + = 0. No assumption 
is made on the ingoing null expansion 9~ of a MOTS. If E is contained in a time 
symmetric Cauchy surface, then 9 + = 0 if and only if E is minimal. 

Marginally trapped surfaces are of central importance in general relativity, where 
they play the role of apparent horizons, or quasilocal black hole boundaries. The 
conjectured Penrose inequality, proved in the Riemannian case by Huisken and 
Ilmanen [HI01] and Bray [BraOl], may be formulated as an inequality relating the 
area of the outermost apparent horizon and the ADM mass. The technique of 
excising the interior of black holes using apparent horizons as excision boundaries 
plays a crucial role in current work in numerical relativity, where much of the 
focus is on modelling binary black hole collisions. 

In spite of the importance of marginally trapped surfaces in the geometry of 
spacetimes, the extent of our knowledge of the regularity and existence of these 
objects is rather limited compared to the situation for minimal surfaces. 

A smooth marginally outer trapped surface is stationary with respect to variations 
of area within its outgoing null cone, in view of the formula 

dfl+^T, = s 

where / is a function on E. The second variation of area at a MOTS in the 
direction l + is 

V» + = -(lx + | 2 + G(i+, (+))/ 

where G denotes the Einstein tensor of L, and x + is the second fundamental 
form of E with respect to l + . For minimal surfaces in a Riemannian manifold, or 
maximal hypersurfaces in a Lorentz manifold the second variation operator is an 
elliptic operator of second order. In contrast, the above equation shows that the 
second variation operator for area of a MOTS, with respect to variations in the 
null direction l + . is an operator of order zero. Therefore, although MOTS can be 
characterized as stationary points of area, this point of view alone is not sufficient 
to yield a useful regularity result. In spite of this, as we shall see below, there is 
a natural generalization of the stability condition for minimal surfaces, as well as 
of the regularity result of Schoen, Simon, and Yau, to marginally outer trapped 
surfaces. 

It is worth remarking at this point that if we consider variations of area of spacclike 
hypersurfaces in a Lorentz manifold, the stationary points are maximal surfaces. 
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Maximal surfaces satisfy a quasilinear non-uniformly elliptic equation closely re¬ 
lated to the minimal surface equation. Dne to the fact that maximal hypersur¬ 
faces are spacelike, they are Lipschitz submanifolds. Moreover, in a spacetime 
satisfying the timelike convergence condition, every maximal surface is stable. 
Hence, the regularity theory for maximal surfaces is of a different flavor than the 
regularity theory for minimal surfaces, cf. [Bar84], 

Assume that L is provided with a reference foliation consisting of spacelike hyper¬ 
surfaces {M t }, and that E is contained in one of the leaves M of this foliation. 
Let ( g , K ) be the induced metric and the second fundamental form of M with 
respect to the future directed timelike normal n. Further, let v be the outward 
pointing normal of E in M and let A be the second fundamental form of E with 
respect to v. After possibly changing normalization, l± — n ± n, we have 

9 ± = H ± tr s K 

where H = tr A is the mean curvature of E and tr^ K is the trace of the projection 
of K to E. Thus the condition for E to be a MOTS, 6 + = 0, is a prescribed mean 
curvature equation. 

The condition that plays the role of stability for MOTS is the stably outermost 
condition, see [AMS05, New87]. Suppose E is contained in a spatial hypersurface 
M. Then E is stably locally outermost in M if there is an outward infinitesimal 
deformation of E, within M which does not decrease 6 + . This condition, which 
is equivalent to the condition that E is stable in case M is time symmetric, 
turns out to be sufficient to apply the technique of [SSY75] to prove a bound 
on the second fundamental form A of E in M. In contrast to the situation for 
minimal surfaces the stability operator defined by the deformation of 6 + is not 
self-adjoint. Nevertheless, it has a real principal eigenvalue with a corresponding 
principal eigenfunction which does not change sign. 

The techniques of [SSY75] were first applied in the context of general relativity 
by Schoen and Yau [SY81], where existence and regularity for Jang’s equation 
were proved. Jang’s equation is an equation for a graph in N = M x R, and is 
of a form closely related to the equation 6 + = 0. Let u be a function on M, and 
let K be the pull-back to N of K along the projection N —> M. Jang’s equation 
is the equation 




DiDjU 

Vl + \Du \ 2 



0 


where g = g lJ — is the induced metric on the graph E of u in N. Thus 

Jang’s equation can be written as 9 = 0 with 


9= H + tr^K, 


where H is the mean curvature of E in N. This shows that Jang’s equation 9 = 0 
is a close analog to the equation 9 + = 0 characterizing a MOTS. Solutions to 
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Jang’s equation satisfy a stability condition closely related to the stably outermost 
condition stated above. This is due to the fact that Jang’s equation is translation 
invariant in the sense that if u solves Jang’s equation, then also u + c is a solution 
where c is a constant. Thus, in the sense of section 5, graphical solutions to 
Jang’s equation are stable. This fact allows Schoen and Yau [SY81] to apply 
the technique of [SSY75] to prove regularity for solutions of Jang’s equation. It 
is worth remarking that although the dominant energy condition is assumed to 
hold throughout [SY81], in fact the proof of the existence and regularity result 
for solutions of Jang’s equation presented in [SY81] can be carried out without 
this assumption. In the present paper, the dominant energy condition is not used 
in the proof of our main regularity result, cf. Theorem 1.2 below. 

It was proved by Galloway and Schoen [GS05], based on an argument for so¬ 
lutions of the Jang’s equation in [SY81], that the stability of MOTS implies a 
“symmetrized” stability condition, which states that the spectrum of a certain 
self-adjoint operator analogous to the second variation operator for minimal sur¬ 
faces is non-negative. The fact that stability in the sense of stably outermost 
implies this symmetrized version of stability was used in [GS05] to give condi¬ 
tions on the Yamabe type of stable marginal surfaces in general dimension. It 
turns out that this weaker symmetrized notion of stability is in fact sufficient 
for the curvature estimates proved here. The symmetrized notion of stability is 
also used in our local area estimates. However, since this notion has no direct 
interpretation in terms of the geometry of the ambient spacetime, we prefer to 
state our results in terms of the stably outermost condition. 

Statement of Results 

The stability condition for MOTS which replaces the stability condition for min¬ 
imal surfaces and which allows one to apply the technique of [SSY75] is the 

following. 

Definition 1.1. E is stably outermost if there is a function f > 0 on E, / ^ 0 
somewhere, such that M + > 0. 

When there is no room for confusion we will refer to a stably outermost MOTS 
simply as a stable MOTS. This is analogous to the stability condition for a min¬ 
imal surface N C M. The condition that there exist a function / on N, f > 0, 
/ 7^ 0 somewhere, such that 5f v H > 0 is equivalent to the condition that N is 
stable. 

The main result of this paper is the following theorem, cf. theorem 6.10, corollary 
6.11 as well as theorem 7.1. 

Theorem 1.2. Suppose E is a stable MOTS in (. M,g,K ). Then the second fun¬ 
damental form A satisfies the inequality 

Plloo < PHA'Iloo, IIVA'IU, || M Rm|| 00 ,inj(M,^)- 1 ). 
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Here || • ||oo denotes the sup-norm of the respective quantity, taken on E. As an 
application we prove a compactness result for MOTS, cf. theorem 8.1. 

Theorem 1.3. Let (g n , K n ) be a sequence of initial data sets on a manifold M. 
Let ( g , K) be another initial dataset on M such that 

|| M Rm|| 00 < C , 

ll^lloo + W m \7K\\ 00 <C, 

inj ( M ,g) > C~ l , 

for some constant C. Assume that 

g n ~> g in Cf oc (M, g) and , 

K n —* K m ClJM,g). 

Furthermore, let E n C M be a sequence of immersed surfaces which are stable 
marginally outer trapped with respect to (g n , K n ) and have an accumulation point 
in M. In addition, assume that the E n have uniformly locally finite area, that is, 
for all x G M there exists 0 < r — r(x) and a = a(x) < oo such that 

|E n 0 Bu tn (x, r)| < ar 2 uniformly in n, 

where BM tn (x,r ) denotes the ball in M around x with radius r. 

Then a subsequence of the E n converges to a smooth immersed surface E locally 
in the sense of C 1,a graphs. E is a MOTS with respect to ( g,K ). If E is compact, 
then it is also stable. □ 

Outline of the paper 

In sections 2 and 3 we discuss the notation and preliminary results, as well as a 
Simons identity which holds for the shear of a MOTS. Section 4 introduces the 
linearization of the operator 9 + acting on surfaces represented as graph over a 
MOTS. The stability conditions we use are discussed in section 5. The curvature 
estimates are derived in section 6 under the assumption of local area bounds. 
In section 7 we show how these bounds can be derived in terms of the ambient 
geometry. Finally section 8 uses the established curvature bounds to prove the 
compactness theorem. 


2 Preliminaries and Notation 


In this section we set up notation and recall some preliminaries from differential 
geometry. In the sequel we will consider two-dimensional spacelike submanifolds 
E of a four dimensional manifold L. As a space time manifold, L will be equipped 
with a metric h of signature (—,+,+,+). The inner product induced by h will 
frequently be denoted by (•, •). In addition, we will assume, that E is contained in 
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a spacelikc hypersurface M in L. The metric on M induced by h will be denoted 
by g, the metric on E by 7 . We will denote the tangent bundles by TL,TM, 
and T E, and the space of smooth tangential vector holds along the respective 
manifolds by T(E), X(M), and X(L). Unless otherwise stated, we will assume 
that all manifolds and fields are smooth. 

We denote by n the future directed unit timelike normal of M in L, which we will 
assume to be a well defined vector held along M. The normal of E in M will be 
denoted by u, which again is assumed to be a well defined vector held along E. 

The two directions n and v span the normal bundle J\f E of E in L , and moreover, 
we can use them to define two canonical null directions, which also span this 
bundle, namely / ± := n ± u. 

In addition to the metrics, h and its Levi-Civita connection L V induce the second 
fundamental form K of M in L. It is the normal part of L V, in the sense that 
for all vector helds X, Y E X(M) 

L VxY = M XxY + K(X, Y)n . (2.1) 

The second fundamental form of E in M will be denoted by A. For vector helds 
X, Ye T(E) we have 

M XxY = "VxY - A(X, Y)u . (2.2) 

For vector helds X, Y E T(E), the connection of L therefore splits according to 

L XxY = s Vv cY + I^(X, Y)n - A(X, Y)v = s Vv :Y - *(X, Y ), (2.3) 

where I(A", Y) = A(X, Y)u — K^(X, Y)n is the second fundamental form of E in 
L. Here K T denotes the restriction of K to T E, the tangential space of E. 

The trace of 1 with respect to 7 , which is a vector in the normal bundle of E, is 
called the mean curvature vector and is denoted by 


H = ^I(e,,e,), (2.4) 

i 

for an orthonormal basis e\, e 2 of E. Since 7i is normal to E, it satisfies 

H = Hu - Pn (2.5) 

where H = ^A^ is the trace of A and P = the trace of A' s , with respect 

to 7 . For completeness, we note that the norms of II and 7i are given by 

|H | 2 = \A\ 2 - |/l s | 2 and (2.6) 

\H\ 2 = H 2 -P 2 . (2.7) 


Recall that since Ti. and II have values normal to E, the norms are taken with 
respect to h and are therefore not necessarily nonnegative. 
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We use the following convention to represent the Riemann curvature tensor s Rm, 
the Ricci tensor s Rc, and the scalar curvature s Sc of E. Here X, Y,U,V G A(E) 
are vector fields. 

s Rm(X, Y, U, V) = ( s Vx S VyC/ - S Vy S VxH - S V[x ,Y]U, V > , 
s Rc(X, Y) = Y S Rm(X, ej, e*, Y ), 

i 

s Sc = 'Y / S Rc(ej, e*). 

i 

Analogous definitions hold for M Rm, M Rc, M Sc and L Rm, L Rc, L Sc, with the ex¬ 
ception that for l Rc and L Sc we take the trace with respect to the indefinite 
metric h. 

We recall the Gauss and Codazzi equations of E in L , which relate the respective 
curvatures. The Riemann curvature tensors s Rm and L Rm of E and L respec¬ 
tively, are related by the Gauss equation. For vector fields X , Y, U, V we have 

s Rm(A, Y, U, V) = 

L Rm(X, Y, U, V ) + (JL(X, V), JL(Y, U)) - (JL(X, U), E(F, V)) . (2.8) 

In two dimensions, all curvature information of E is contained in its scalar cur¬ 
vature, which we will denote by s Sc. The scalar curvature of L will be denoted 
by L Sc. The information of the Gauss equation above is fully contained in the 
following equation, which emerges from the above one by first taking the trace 
with respect to Y, U and then with respect to X, V 

s Sc = L Sc + 2 L Rc(n, n ) — 2 L Rc(z/, v) — 2 L Rm(iy n, n, u) + \H\ 2 — |II| 2 . (2.9) 

The Codazzi equation, which relates L Rm to E, has the following form 

< L VxE(R, Z),S) = < V y E(X, Z), S) + L Rm(X, Y, S, Z) (2.10) 

for vector holds X , Y, Z G A(E) and S G T(A/’E). 

There is also a version of the Gauss and Codazzi equations for the embedding of 
M in L. They relate the curvature L Rm of L to the curvature M Rm of M. For 
vector fields X, Y, U, V G X(M) we have 

M Rm (X,Y,U,V) 

= L Rm(X, Y, U, V ) - K(Y, U)K(X, V ) + K(X, U)K(Y, V ), 1 J 

m V x K(Y, U) - M XyK{X , U) = L Rm(X, 1", n, U). (2.12) 

These equations also have a traced form, namely 

M Sc = L Sc + 2 l Rc (n,n) — (tr K) 2 + \K| 2 and (2-13) 

M div/l - M V tr K = L Rc(-, n) . (2.14) 
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We now investigate the connection A V on the normal bundle WE of E. Recall 
that for sections N of A/"E and X e T(E), this connection is defined as follows 


N XxN = ( i VxiV) ± , 

where again (•)-*- means taking the normal part. We have 
0 = X(l) = X((n,n}) = 2( N V x n,n ), 
and similarly ( w Vxg v) = 0. Therefore the relevant component of A V is 
< A Vxgn) = ( l Va :v,n) = -K(X, v ). 

Recall that X is tangential to E. This lead us to define the 1-form S along E by 
the restriction of K(-, v) to T E. 

S(X):=K(X,u). (2.15) 

Then, for an arbitrary section N of WE with IV = fv + gn, we have 
N XxN = X{f)u + X(g)n + S{X)(fn + gv ). 

In particular 

N Xxl ± = ±S(X)l ± . (2.16) 

We will later consider the decomposition of II into its null components. For 
X, Y e T(E) let 

X ± (X, Y) := (I{X, Y), l±) = K(X, Y) ± A(X, Y). (2.17) 

The traces of y 1 * 1 respectively will be called 6 ± 

0±= {H,l ± ) = P±H. (2.18) 

The Codazzi-equation (2.10) implies a Codazzi equation for y ± . 

Lemma 2.1. For vector fields X, Y, Z e T(E) the following relation holds 

V xX ±(Y,Z) = Vyx ± (X,Z)+Q ± (X,Y,Z)tx ± (X,Z)S(Y)±x ± (Y,Z)S(X) . 

(2.19) 

Here 


Q ± (X, Y, Z ) = L Rm(X, Y, l ± ,Z). 


( 2 . 20 ) 



3 A Simons identity for % ± 


We use the Codazzi equation we derived in the previous section to compute an 
identity for the Laplacian of y ± , which is very similar to the Simons identity for 
the second fundamental form of a hypersurface [Sim 68 , SSY75]. 

The Laplacian on the surface E is defined as the operator 
S A = • 

In the sequel, we will drop the superscript on S A and S V, since all tensors below 
will be defined only along E. We will switch to index notation, since this is 
convenient for the computations to follow. In this notation 

1 h-3 q 

denotes a (p, g)-tensor T as the collection of its components in an arbitrary basis 
for the tangent spaces. To make the subsequent computations easier, we 
will usually pick a basis of normal coordinate vectors. Also note that we use Latin 
indices ranging from 1 to 2 to denote components tangential to the surface E. 

Recall, that the commutator of the connection is given by the Riemann curvature 
tensor, such that for a ( 0 , 2 )-tensor 

— VzVfcTy = S R ITiklmiTmj + S Rm klmjTim • (3.1) 

Note that we use the shorthand s R.rn fcimj T, m = T, Rm k i p jT iq 'y pq , when there is 
no ambiguity. That is, we assume that we are in normal coordinates where 
7 ij = y *- 7 = Sij. Also note that this fixes the sign convention for E R,m. y fc; such 
that s Rcjj = ^Ruiikkj is positive on the round sphere. 

Lemma 3.1. The Laplacian of x — X + satisfies the following identity 

Xij^Xij Xq T Xij ( R ^LVkHkXlj Y R X&kiljXkl) 

Y Xij^ki^Qkij XkjSi T Xij^k) T Xij^i^Qkjk @ Sj T XjkSk) 

- |i| 2 lxl 2 + 0+xtjXtkXti - o + xtjXp<l - PxtjX+kXti 

where P = ^ Kfj is the trace of /X s . 

Proof. Recall that in coordinates the Codazzi equation (2.19) for Xij reads 

^ iXjk V jXik T Qijk Xik^j T Xjk^i ■ (3.2) 

Then compute, using (3.2) in the first and third step, and the commutator relation 
(3.1) in the second, to obtain 

VfcV iXij ^ k^iXlj Y ^ki^Qlij Xij Y Xij^l) 

VjVfcXp Y bill kjjiil\arj A RaikimjXlm 

Y Vk(Quj — Xij^i + XijSi) (3-3) 

Vj VjXA:/ Y Rm kimlXmj Y Rm kimjXlm 

Y V ki^Qlij XijSi Y Xij^l) Y V i(^Qkjl Xkl^j Y Xjl^k) ■ 
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We will use the Gauss equation (2.8) to replace the s Rm-terms by ^Rm-terms. 
Observe, that 


n o = ~\x% 1 ~ \xj 


2 A-ij 


Plugging this into the Gauss equation (2.8) gives 

SR m ijki = L Rm ijk i + ^(xtkXji + XikXji ~ xtiXjk ~ XuX %) • 

Combining with (3.3), we infer that 


- rv c i i ov /v.' 


V/.Y/x,, = V,Y;\/,/ • 

+ HXiiXfem + XaXkm ~ XklXim 
+ KXfcmXij + XkmXtj ~ XkjXim 
~t~ Y ki^Qlij 


Vmj 


XklXtm) Xr, 

j ~~l"~ X'kmX'ij X'kjX'im - x kj x im )x] m 
Xlj^i T T Y i^Qkjl XklSj + Xjl^k) 


Taking the trace with respect to k, l yields 

A Xij Y/ Yy T bill kllkXIj T hill k/lj X hi 

Y k (Q /.■?/ XkjSi T Xij^k) T Y ii^Qkjk @ Sj T XjkSk ) 

+ ^(x£|x + | 2 + (x + , X _ )xS - 0 + X%Xki ~ o~x%xi) 

+ \{x%Xkixti - xJfcXfcixii) 

We contract this equation with yh and obtain 

Xij^Xij A/./ Y/ Yy 0 T Xij ( bui k ! /ikXij T RmfctZjXfc/) 

T Xij^ki^Qkij XkjSi T Xij^k) T XijY ii^Qkjk @ Sj T XjkSk ) 

+ (x + , x')lx| 2 - ¥ + xtjX%Xki - ¥~xijX%Xu ■ 

Now observe that = 2AA — yA and = 2P — 0 + . Substituting this into 
the last two terms, together with (y + ,x ) = — |II| J , we arrive at the identity we 
claimed. □ 


4 The Deformation of 9 + 


This section is concerned with the deformation of the operator 6 + , as defined 
in equation (2.18). We begin by considering an arbitrary, spacclikc hypersurface 
S C L. Assume that the normal bundle is spanned by the globally defined null 
vector fields l ± , such that (/ + , l~) = —2. We call such a frame a normalized null 
frame. As before, let 9 ± := (H, l ± ). We abbreviate x — X + - 

A variation of S is a differentiable map 

F : S x (—£, e) —> L : (x, t) t—■> F(x, t ), 
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such that F(-, 0) = ids is the identity on E. The vector field = V is called 

variation vector field of F. We will only consider variations, with variation vector 
fields V of the form V = al + + (3l~. 

Note that in this setting, as a normalized null frame is not uniquely defined by 
its properties, the notion of 9 + depends on the frame chosen. The freedom we 
have here is the following. Assume AA is another normalized null frame for the 
normal bundle of E, that is Ai(A; , AA) = 0 and h(k + ,k~) = —2. Since the null 
cone at each point is unique, the directions of AA can be aligned with l ± . But 
their magnitudes can be different, so k + = e u l + and k~ = e~ u l~ with a function 
u e C°°( E). 

Therefore, if we want to compute the deformation of 9 + , it will not only depend 
on the deformation of E, as encoded in the deformation vector V. It will also 
depend on the change of the frame, that is on the change of the vector l + , which 
is an additional degree of freedom. 

To expose the nature of that freedom, observe that if P(A) is a null frame on each 
E t := F(E,t ), then ^ is still normal to E. On the other hand 

0= SI, = „{i + .' + ) = 2(ai| i _ o .a> and 
°= sl<.o< (+ ’ r > = < ^r| t , 0 ' r > + ( fr|,_ 0 - (+ ) 

Therefore ^ = wl± for a function w € C'°°(E). Thus the linearized change 

of the frame is described by the single function w, which we will call the variation 
of the null frame. 

If we fix both of the quantities V and w, a straight forward (but lengthy) com¬ 
putation gives the deformation of 6 + . 

Lemma 4.1. Assume F : Ex (—e,e) —> L is a variation of E with variation 
vector field V = al + + (3l~. Assume further that the variation of the null frame 
is w. Then the variation of 9 + is given by 

S V}W 9 + = 2A (3 - AS{V(d) - a(\ X \ 2 + L Rc(/+/+)) + 2 9 + w 

— /3(2divS' — 2|S'| 2 — |E| 2 + L Rc(/+/-) - | L Rm {Ifl~, l~,l + )). 

If we consider marginally trapped surfaces, then the term 9 + w in the previous 
calculation vanishes, and we get expressions independent of the change in the 
frame. As a consequence, we state the following two corollaries, which also restrict 
the variations we take into account. 

Corollary 4.2. Assume E is a marginally trapped surface, that is, it satisfies the 
equation 9 + = 0. Then the deformation of 9 + in direction of —l~ is given by 

fi-p l -,v>9 + — 2L_/3, 
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where the operator A_ is given by 


L_(5 = -A (3 + 2 A(V/3) + /3(div S - ||H| 2 - |A | 2 - _), 
and = | L Rm(/+/-/-/+) - | L Rc(/+/-). 

If we assume that E C M, where M is a three dimensional spacclikc surface, then 
E can be deformed in the direction of u, the normal of E in M. The deformation 
of 9 + then turns out to be the following. 

Corollary 4.3. Assume E is a marginally trapped surface, then the deformation 
of 9 + in the spatial direction of v := |(/ + — l~) is given by 

^ fu,w -Lm f > 

where the operator L M is given by 

L u f = -A/ + 2A(V/) + /(div A - | X | 2 + <A S , x ) - |5 | 2 - T M ) , 
and d> M = | L Rm (lfl~l~l + ) + L Rc(u,l + ). 

Remark 4.4. (i) Using the Gauss equation (2.9), we can rewrite the expression 
for L m as follows 

L M f = -A/ + 2S(V/) + /(div5 - i | X | 2 - |5 | 2 + f Sc - T M ) • (4.1) 

Here A! M — G(n,l + ) where G = L Rc — | L Sc h denotes the Einstein tensor of h. 
Note that in view of the Gauss and Codazzi equations of the embedding M ^ L , 
equations (2.13) and (2.14), the term can be rewritten as 

= i( M Sc + (trA) 2 -|Af)-( M divA'- M VtrA» =: 8 tt(//- J{u)) , (4.2) 

where 8nJ = M divK — A/ V tr K is the projection of G(n,-) to M and 167 t/x = 
M Sc + (tr A ') 2 — |A '| 2 = G(n,n). The dominant energy condition is equivalent 
to \J\ < fi. Thus, if the dominant energy condition holds, turns out to be 
non-negative. 

(ii) The same procedure gives that we can write L_ as 

L-f = -A/ + 2S(Vf) + / (div A - |A | 2 + fSc - T_) . (4.3) 

with = G(l + ,l~). Note that is non-negative if the dominant energy 
condition holds. However, this representation does not contain a term |x| 2 , which 
does not allow us to get estimates on sup |x| 2 - However, in the case of strict L_ 
stability there is a sheet M such that the surface is A^-stable. We can then apply 
the subsequent results to get the estimates of theorem 1.2 in this case. □ 
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5 Stability of marginally outer trapped surfaces 


As before, consider a four dimensional space time L 4 , with a three dimensional 
spacelikc slice M 3 . As in the previous sections, the future directed unit normal 
to M in L will be denoted by n. In M consider a two dimensional surface E, 
such that there exists a global unit normal vector field v of E in M. The vector 
fields n and v span the normal bundle of E in L and give rise to two canonical 
null vectors Z* = n ± v. Again we use the shorthand x = X + ■ 

In this section we will introduce two notions of stability for a marginally trapped 
surface. These are related to variations of the surface in different directions. The 
first definition is equivalent to definition 2 in [AMS05]. There a stably outermost 
marginally outer trapped surface, is defined as surface, on which the principal 
eigenvalue of Lm is positive. Here an Lm- stable MOTS is defined as follows. 

Definition 5.1. A two dimensional surface E C M C L is called a L M -stable 
marginally outer trapped surface if 

(i) E is marginally trapped with respect to l + , that is 9 + = 0. 

(ii) There exists a function f > 0, / ^ 0 such that Lm f > 0 . Here Lm is the 
operator from corollary f.3. 

Remark 5.2. (i) Although Lm is not formally self-adjoint, the eigenvalue of Lm 
with the smallest real part is real and non-negative (cf. [AMS05, Lemma 1]). 
This definition is equivalent to saying, that the principal eigenvalue of Lm is 
nonnegative. This is seen as follows: 

Let A be the principal eigenvalue Lm- Then, since A is real, the L 2 -adjoint L* M of 
Lm has the same principal eigenvalue and a corresponding eigenfunction g > 0. 
Pick / > 0 as in the definition of Z/M-stability, ie. Lm/ > 0. Then compute 

A / fg dp= [ fL* M g dn= [ L M fg d/i . 

,/E JT, J S 

As / > 0, / ^ 0, g > 0 and L M f > 0, this implies A > 0. 

The eigenfunction if of L M with respect to the principal eigenvalue does not 
change sign. Therefore it can be chosen positive, if > 0. Thus, the definition 
in fact is equivalent to the existence of if > 0 such that L M if = Xif > 0. We 
will use this fact frequently in the subsequent sections. Note that LM-stability is 
equivalent to the notion of a stably outermost MOTS in [AMS05, Definition 2]. 
(ii) The conditions from the above definition are satisfied in the following sit¬ 
uation. Let E = dLl be the boundary of the domain Z2 and satisfy 6 + = 0. 
Furthermore assume that there is a neighborhood U of E such that the exterior 
part U\il does not contain any trapped surface, ie. a surface with 9 + < 0. Then 
E is stable. Assume not. Then the principal eigenvalue would be negative and 
the corresponding eigenfunction if would satisfy Lm/ < 0, if > 0. This would 
imply the existence of trapped surfaces outside of E, since the variation of E in 
direction ifv would decrease 9 + . □ 
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Note that the condition 6 + = 0 does not depend on the choice of the particular 
frame. Therefore, to say that a surface is marginally trapped, we do not need 
any additional information. In contrast the notion of stability required here does 
depend on the frame, since clearly there is no distinct selection of v when only E 
— and not M — is specified. 

To address this issue, we introduce the second notion of stability of marginally 
outer trapped surfaces, namely with reference to the direction —l~. This definition 
is more in spirit of Newman [New87] and recent interest in the so called dynamical 
horizons [AK03, AG05]. 

Definition 5.3. A two dimensional surface E C M C L is called a L_-stable 
marginal outer trapped surface (L - -stable MOTS) if 

(i) E is marginally trapped with respect to l + , that is 9 + = 0. 

(ii) There exists a function f > 0, / ^ 0 such that L~f > 0. Here is the 
operator from corollary f.2. 

Remark 5.4. It turns out that this notion of stability does not depend on the 
choice of the null frame. This is due to the natural transformation law of the 
stability operator L_ when changing the frame according to l + = fl + and l~ = 
/ _1 / _ . Then the operator with respect to this frame satisfies f~ 1 L(f/3 ) = LB 
for all functions [3 G C'°°(E), as it is expected from the facts that 6 + = f6 + and 
-pi- = -pfl~. □ 

Remark 5.5. (i) Remark 5.2 is also valid here, in particular the definition implies 
that there exists a function if > t) with L_if > 0. 

(ii) Technically speaking, the equation for a marginally trapped surface pre¬ 
scribes the mean curvature H of E in M to equal minus the value of a function 
P : TM —> R : (p, v) i—> trK — K l3 i/v\ namely H{p) = —P(p, v) for all p G E. 
This is a degenerate quasilincar elliptic equation for the position of the surface. 
These equations do not allow estimates for second derivatives without any ad¬ 
ditional information. This is where the two stability conditions come into play. 
They give the additional piece of information needed in the estimates as in the 
case for stable minimal surfaces. □ 

The two notions of stability above imply the positivity of certain symmetric dif¬ 
ferential operators as it was noticed in [GS05] for the operator L M . However, 
the inequality there is not quite sufficient for our purposes, it needs some further 
rearrangement. This is the content of the following Lemmas. Basically lemmas 
5.6 and 5.6 are the only way how stability is used in the subsequent estimates. 
Basically one could use these, in paricular equation (5.1) to define a notion of 
symmetrized stability for MOTS. 

Lemma 5.6. If E is a stable MOTS , then for all £ > 0 for all rj G C'^°(E) the 
following inequality holds 

[ V 2 \x\ 2 dp < (1 +e) / |V//| 2 + r/ 2 ((4e)" 1 |A' E | J ~ ^m) dpi. 

J E J£, 
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Proof. Take / as in the definition of a stable MOTS. From remark 5.2 we can 
assume / > 0. Then f~ l L M f > 0. Multiply this relation by rf, integrate, and 
expand L M as in corollary 4.3. This yields 

0 < f V 2 ( - / _1 A/ + 2/- 1 S(V/) + div S - \x\ 2 + (A s , X ) ~ l^| 2 - v I ; m) d/i. 

By sorting terms, and partially integrating the Laplacian and the divergence term, 
we obtain 

J yixi 2 + y(/- 2 iv/r - 2 /- i s(v/) + isi 2 ) dj, 

< f .r'V/ - S> + y Ixl | A' E | - i/ 2 'I'm d/i . 

By the Schwarz inequality 

[ 2r/(V/ 7 , / _1 V/ — S') d/i < [ \Vitf + V 2 \f- l Vf - S\ 2 d/i 
and for any £ > 0 

/‘r/ 2 |iF s || X |d/i<(4 £ )- 1 f r/ 2 |iF s | 2 d/i + £ /Vlx| 2 - 

Cancelling the terms / s ? 7 2 |/ _1 V/ — S'pd/i and £ J s -// 2 |x| 2 d/i on both sides, we 
conclude the claimed inequality. □ 

The following lemma is based on the original computation of [GS05]. 

Lemma 5.7. If £ is a stable MOTS, then for all e > 0, there exists C(e _1 ) such 
that for all rj G C£°(£) the following inequality holds: 

f r] 2 \x\ 2 dfi < (1 + e) I |V//| 2 + t/ 2 (| m Rc| - + C(e~ l )\K s \ 2 ) d/i. 

Jt. J s 

Proof. We proceed as in the computation of lemma 5.6, but with the alternative 
representation (4.1) for L M . As in [GS05], we get 

[ 77 2 |x| 2 d/i < [ 2|V//| 2 + r/ 2 ( s Sc - 2T m ) d/x. (5.1) 

We use Gauss’ equation to replace E Sc in the following way 
E Sc = M Sc - 2 m Rc(u, u) + H 2 - \A\ 2 , 

where A is the second fundamental form of £ C M and H is the mean curvature. 
We can move the |A| 2 term to the left hand side. Then H 2 = P 2 by 6 + = 0 and 
thus H 2 < 2|A' S | 2 . The remaining terms are controlled by j u Rc|. Inserting this, 
we find that 

[ V 2 (\x\ 2 + |A| 2 ) d/i < 2 f |Vr/| 2 + r/ 2 (| M Rc| + |AT S | 2 - tf M ) d/i. 

J s ,/s 
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Now fix £ > 0. Since x = -A + Zl s we can estimate 

2|*| 2 < (1 + £)W 2 + (1 - e)(|A| 2 + 2{A, K*) + \K*\ 2 ) 

< (1 + e)(\x\ 2 + \A\ 2 ) + (2 - 2e + (2e) _1 )|A' s | 2 . 

Inserting this into the above inequality we find the claimed inequality. □ 

A similar, but fundamentally different inequality holds in the case of Testability. 
The fundamental difference is that the gradient term on the right has a factor of 
a little more than two, instead of a little more than one, as with the Z/M-stability. 
In view of lemma 5.9 this factor of two is not at all surprising. This factor is the 
reason why the procedure in section 6 does not work to give curvature estimates 
for L_-stable surfaces. 

Lemma 5.8. If E is a L--stable MOTS, then for alls > 0 the following inequality 
holds: 

[ IxlVd// < 2(1 + e ) [ | V 77| 2 + 77 2 (( 2 £) _ 1 | A : s | 2 - d/x. 

JT. J S 

We conclude with the remark that Z/M-stability implies Testability. 

Lemma 5.9. Let ( L,h ) satisfy the null energy condition, i.e. assume that for all 
null vectors l we have that L Rc(/,/) > 0. Then if E is an L M -stable MOTS, then 
it is also L--stable. 

Proof. We use the notation from section 4, where we introduced the linearization 
of 9 + . For any function / compute 

L M f - L-f = 5 fu , w 6 + - \8fi-, w 0 + = \5 fl+ , w 0 + = -\f{\ X \ 2 + l Rc(/ + , 1+)). 

If / > 0, then by the null energy condition, the right hand side is non-positive. If 
in addition Lm/ A 0, as in the definition of Z/M-stability, then this implies that 

L-f > L M f > 0. 

Hence E is also Lm stable. □ 


6 A priori estimates 

In this section we derive the actual estimates for stable outermost marginally 
trapped surfaces. All but the most basic estimates hold only for LM-stable sur¬ 
faces, as defined in section 5. This is due to the factor of two appearing in front 
of the gradient term in lemma 5.8, which does not allow us to carefully balance 
the Simons inequality and the stability inequalities. 

Throughout this section we will make the assumption, that the surfaces in ques¬ 
tion have locally uniformly finite area. 
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Definition 6.1. If there exists r > 0 and a < oo, such that for all x G £ 

\Bz(x,r) \ < a, (6.1) 

then we say that £ has (r, a)-locally uniformly finite area. 

Here B^(x, r) denotes the the ball of radius r around x in £. In the sequel we 
will denote B^(x, r) by B(x,r). The estimates below work in exactly the same 
way if intrinsic balls are replaced by extrinsic balls. Later, we will derive such 
bounds for stable MOTS in terms of the ambient geometry. 

We first begin with the observation that the stability of a MOTS gives a local 
L 2 -estimate for the shear tensor x — X + - 

In the sequel, for a tensor T, we denote by HTHoo = sup E \T\. That is, oo-norms 
are taken on £ only. Constant are always denoted by C, if we want to clarify the 
dependence of the constants, we denote by C(a, b,.. .) a constant that depends 
on the quantities a, b... in such a way that C deteriorates as a + b + oo. 

Lemma 6.2. Suppose £ is an Lu-stable MOTS with (r, a)-locally uniformly finite 
area. Then for all x G £ 

I \x \ 2 d/i < C(r~%a, ||iC s ||oo, ||'b A r||oo) • 

J B{x,r/ 2) 

Alternatively, the constant can be chosen to depend on ||^m||oo an d || M Rc|| 00 
instead of ||'L M || 00 . 

Proof. The desired bound is easily derived from lemma 5.6 or 5.7. We will restrict 
ourselves to the proof of the first statement. To this end, fix £ — |, x e £ and 
choose a cut-off function rj > 0 such that rj = 1 on B(x,r/ 2), rj — 0 on dB(x,r), 
and |V//j < 4r _1 . The left hand side of the equation in lemma 5.6 is then is an 
upper bound for the left hand side in the claim, whereas the right hand side can 
be estimated by the claimed quantities. □ 

This estimate can also be derived from L_-stability: 

Lemma 6.3. Suppose £ is an L^-stable MOTS with (r, a)-locally uniformly finite 
area. Then for all x G £ 


\x\ 2 dp < C(r \a, H-R^Hoo, ||'L_|| 0 o) • 


J B(x,r/ 2) 

Proposition 6.4. Let £ be an Lm- stable MOTS. For any e > 0, any p > 2, and 
any function rj we have the estimate 

[ V 2 \x\ P+2 dp < ^(l + £) [ V 2 \x\ P ~ 2 \V\x\\ 2 dji 
J s ds 

+ C{e~\ Il'MU Halloo) J (- n 2 + |V7 7 | 2 ) \x\ p dji. 

Alternatively, we can make the constant on the right hand side to be of the form 
C'(£- 1 ,||T m || 00 ,||A^|| 00 ,|| m Rc|| 00 )). 
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This proposition also holds for L_-stable surfaces with appropriate modifications 
of the dependencies of C, and a general factor of two on the right hand side. This 
factor of is the reason, why an argument like the subsequent one fails to give 
curvature estimates for L_-stable MOTS. 

Proof. We will restrict to the proof of the first statement, since the other is proved 
similarly. From lemma 5.6 we find that for any 6 > 0 there is (7(<5 -1 ) such that 
for all functions we have 

[ 0 2 |y| 2 d^i < (1 + 5) I \\/f\ 2 + f 2 (C(5- 1 )\K^\ 2 -^ M )d f i. 

Jt. Jy. 

We substitute f by q\x\ p ^ 2 - To this end compute 

V(/ 7 |xl p/2 ) = V? 7 |y| p/2 + V f |x| p/2_1 V|x|. 


For any 5 > 0 we thus can estimate 

|v(i)|x|" /2 )| 2 < (1 - S) d| x p.-y|V|*|| 2 + c(i- 1 )|v„Hxr. 

Inserting this estimate into the hrst inequality, we find 

[ ri 2 \x\ r+2 <l»< (l+ Sf£ f y|xr 2 |V|x || 2 

J E J S 

+ J C(S~ 1 )\'Vr 1 \ 2 + 7 7 2 | X | p (C'(r 1 )|iF s | 2 - vk M ) d/ 1 . 

Adjusting 5 yields the claimed estimate. □ 

We now aim for an estimate on the gradient term on the right hand side of the 
estimate in proposition 6.4. The main tool will be the Simons identity from 
section 3. To avoid that the estimated depend on derivatives of curvature, we use 
similar techniques as in [Met04], 

Proposition 6.5. Let E be an L M -stable MOTS. Then there exists p 0 > 2 such 
that for 2 < p < p 0 and all functions q we have the estimate 

/yixrMxiiw 

< C(p, ||f M |U ||V E |U, IIQIU, || 1 Rm E |U, ||S|U) /(V + |V„| 2 )|xr + I) 2 |xl p - 2 - 

Alternatively, as before, we can replace the dependence of the constant on ||Tm||oo 
by a dependence on ||Tm||oo || M Rc||oo. 

Before we can start the proof of the proposition, we state the following lemma. 
It states an improved Kato’s inequality similar to [SY81]. A general reference for 
such inequalities is [CGHOO]. 
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Lemma 6.6. On a surface E with 9 + = 0 we have the estimate 

Vxf - |V|xli 2 > i(|V| X || 2 + |Vx| 2 ) -c(\Q\ 2 + |S| 2 |x| 2 ) . 

Here c is a purely numerical constant. 


Proof. The proof goes along the lines of a similar argument in Schoen and Yau 
in [SY81, p. 237], but for the sake of completeness, we include a sketch of it here. 

In the following computation we do not use the Einstein summation convention 

i 12 

and work in a local orthonormal frame for T E. Let T := |Vy| 2 — V|x| • We 
compute 

\x\ 2 T = |x| 2 |Vx| 2 - i|V|x| 2 | 2 

= y ; ixij^kxim.) _ y " ^ y ^ xa^ kx.ij'j 

i,j,k,l,m k ij 

2 ^ ^ kXlm kXij^j 

i,j,k,l,m 

In the last term consider only summands with i = k and j = m. This gives 

\x\ 2 T (xij^iXji - XjNiXij ) > | XijViXii - XjiViXij ) • 

i,j,l l i,j 


Lise the Codazzi equation (3.2) to swap indices in the gradient terms. We arrive 
at 

I x\ 2 T > | (XijViXij+XijQuj-XijQiji)+^2 {QSiXii-Xii'y i 6)-\x\ 2 Si^ . 

I i,j i 


By the fact that (a — b ) 2 > |a 2 — b 2 , this implies 

\x\ 2 T > j 

l i,j 

- I (XijQnj - xijQiji) + xu(0Si - v<0) - |xl 2 ^) 

l i,j i 

>^lxl 2 |v|xl| 2 -c|xr(M 2 + |5| 2 |x| 2 ). 

Dividing by |x| 2 , we get 


IVxl 2 - |V|xlf > isMxll 2 - c(|Q| 2 + |S| 2 !x| 2 ) ■ 


Adding yj(|Vy| 2 — |V|x|| 2 ) to both sides of this inequality and multiplying by 
|| yields the desired estimate. □ 
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Proof of proposition 6.5. We will restrict to the proof of the first statement, since 
the other is proved similarly. Compute 

A|x| 2 — 2|x|A|x| + 2 1V|x| | 2 . 

On the other hand 

A|x| 2 = - + 2|Vx| 2 ■ 

Subtracting these equations yields 

IXI A | x | = Xij^Xij + | Vy| 2 — | V | x 11 2 • 

In the case 9 + = 0, the Simons identity from lemma 3.1 gives 

Xij^Xij Xij ( ^P^-kilkXlj T Km kiljXkl) |E| |x| ^XijXjkXki 

T ki^Qkij XkjSi T Xij^k) T Xij^ii^Qkjk T XjkSk ) • 

Note that XijXjkXki — tr(x 3 ), and the trace of a 2 x 2 matrix A satisfies the 
relation tr A 3 = trA(trA 2 — det A). Since x is traceless, this term vanishes. In 
addition |E| 2 = (x,X~) = lx| 2 - 2(/l s ,x). 

As we are not interested in the particular form of some terms, to simplify notation, 
we introduce the ^-notation. For two tensors Tj and T 2 , the expression Tj * T 2 
denotes linear combinations of contractions of T\ (g) Tj. 

To remember that in the above equation we need to evaluate L Rm only on vectors 
tangential to E, we use the projection of L Rm to TE and denote this by L Rm s . 
Then the above equations combine to 

-IxIAlxl + IVxlHVIxll 2 = |x| 1 + |x| 2 »X*X E +x*X» I Rm s + x»V(Q+x*S) . 

(6.2) 


Multiply this equation by ?/ 2 |x| p 2 and integrate. This yields 

J C — r ? 2 |x| P_1 A|x| +r7 2 |xr 2 (|V X | 2 - |V|x|| 2 ) d/i 

= f h 2 (lxl /,+2 + \x\ p x * A' s + |x| p_2 X * X * L Rm s + |x| p ~ 2 X * V(Q + x * S)) d/i. 

Next, do a partial integration on the term including the Laplacian and on the 
last term on the second line. We find that 

I V 2 (p- l)lxM Vlxl I 2 + ylxr 2 (lVxl 2 - IVlxl I 2 ) 

< [ I) 2 |xl p+2 d(. + c f „\x\ p -'\Vr,\ IVIxll + i)|xr‘|Vi)|(|<?| + |xl|S|)d M 

Jt, 

+ c jjf{\xr l \ K "-\ + lxP| L Rm s | + |xr 2 (|Vx| + |V|xl|)(|Q| + |xl|S|)}d^. 
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(6.3) 


Here c is a purely numerical constant. For any £ > 0, we can estimate 

c f rj 2 \x\ p+1 \K E \dn < £ f p 2 \x\ p+2 dp+ C(£~ 1 ) f r] 2 \x\ p \K^\' 2 dfi 
J s S J s 

as well as 

c/Vlxr 2 (|Vxl + |v|xl|)(IO| + M|S|)d*i 

<£ f r7 2 |xr _2 (| Vx | 2 + | V|xl | 2 ) d p + C^ 1 ) [ r/ 2 (\x\ p \S\ 2 + |y| p_ 2 |<3| 2 ) d./x. 

J £ </£ 

In addition we estimate 

cf ^?lxl p—1 1 Vr^l |v|x|| dn< £ /u 2 |V|xl| 2 | xr 2 + c(e~ l ) I \V V \ 2 \x\ r dn, 

J S «/ S «/ s 


and 


[ oixr'iv^itiei + mis|)<u < [ \vr,\ 2 \ x \* +» 2 (ixr 2 iei 2 +ixnsf) ^ 

7 s J s 

Inserting these estimates into the estimate (6.3), gives 

/>) 2 (p-l)lxr 2 |V|xlf + i? 2 |xr 2 (|Vx| 2 -|V|x|| 2 )df< 

< [ (1 + £)p 2 \x\ p+2 + 2^77 2 1x| p 2 (I Vx | 2 + | V|x11 2 ) d/x + C(e _1 ) / | Vp\ 2 \x\ p d/x 

J s J s 

+ C( £ “M|V !: ||oo, l| L Rm s ||co, ||S||co, IIQIU) J rf(\x\ p + Ixr 2 )dp. 

(6.4) 

We apply lemma 6.6 to estimate the second term on the left hand side from below 
by y; /^(Mxll + |Vx| 2 ) d/x. In addition, use proposition 6.4 to estimate the 
first term on the right hand side. This yields 

^V(p- 1 )lxr _2 | V|x|f + (35 - 2£)r/ 2 |x| p_ 2 (|Vx | 2 + |V|y|| 2 ) d/x 

< fd + e) 2 / s , 2 |xr 2 |V|x|| 2 d, 

+ C{e~\ ||^ S ||oo, IIQIloo, || i Rm s || 00 , Moo) J ri 2 (\x\ p + \x\ P ~ 2 ) + |Vr/| 2 |yf d/x . 

Choose po > 2 close enough to 2 and £ small enough, such that for 2 < p < p 0 
the gradient term on the right hand side can be absorbed on the left hand side. 
This gives the desired estimate. □ 
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Combining propositions 6.4 and 6.5 with the initial Id-estimate in lemma 6.2 
gives the following L p estimates for |y|. 

Theorem 6.7. There exists po > 2 such that for all 2 < p < p 0 and all LiM-stable 
MOTS E which have (r, a)-locally uniformly finite area, the shear x satisfies for 
all x G E the estimates 


[ lxf +2 d^<C. 

' B(x,r/8) 

(6.5) 

f ixr 2 |vw| 2 dj.<c, 

' B(x,r/ 8) 

(6.6) 


and 


/ |V X | 2 d P<C. (6.7) 

J B(x,r/8) 

Here C = (^(r -1 , a, ||^m||oo, Halloo, ||Q||oo, I^Rm^H^, ||S'|| 00 ), or alternatively, 
we can replace the dependence on ||^m||oo by ||'1 / m||oo an d H^RcHoq. 

Proof. Fix x G E and choose rj to be a cut-off function with rj = 1 in B(x,r/ 4) 
and rj = 0 outside of B(x,r/2), such that 0 < rj < 1 and | V?;| < 4r _1 . Plugging 
this into the estimates of 6.5 and 6.4 for p — 2, in view of the local area bound, 
and the local Id-estimate from lemma 6 . 2 , yields Id-estimates for |y| in B(x, r/4). 

We pick po a little smaller, than the value allowed by Proposition 6.5. Then for 
any 2 < p < p 0 , proceed as before, but now choose a cutoff function f) with 
fj = 1 in B(x,r/ 8) and f) = 0 outside of Bfxpr/4), such that 0 < fj < 1 and 
|Vf,| < 8 r 1 . The resulting L p and L p 2 -norms of % on the right hand side 
can now be estimated by combinations of the Id-norm of \x\ and the local area 
bound. 

To see the last estimate, note that in the proof of proposition 6.5, by appropriately 
choosing £, we can retain a small portion of the term f y p 2 |x| p ~ 2 | Vy | 2 d/i on the 
right hand side. □ 

For the next step - the derivation of sup-bounds on y - we use the Hoffman- 
Spruck Sobolev inequality in the following form [HS74], 

Lemma 6.8. For ( M,g) exist constants Cq , cf, such that for all hypersurfaces 
E C M and all functions f G (d 00 (E) with |supp/| < Cq the following estimate 
holds: 

(7l/| 2 dd 7 <cf^|V/| + |/ff|d,i. 

Here H is the mean curvature of E and the constants Cg, cf depend only on a 
lower bound for the injectivity radius and an upper bound for the curvature of 

C M,g ). 
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Remark 6.9. Replacing / by f p in the above inequality and using Holders in¬ 
equality gives that for all 1 < p < oo and all / with |supp/| < Cq 

S CpIsupp/I^^IV/l 2 + |ff/| 2 dp. 

The constant only depends on cf and p. 

Theorem 6.10. Let E be an L M -stable MOTS with (r, a)-locally finite area. Then 
the shear x satisfies the estimate 

sup |x | < C. 
s 

The constant C depends only on r^ 1 , a, ||\I/m||oo; ||/L s ||oo> IIQIIoo? || L R mS || 00; 
||-5||oo, || M Rm|| 0O; and inj(M, g)~ l . 

Alternatively, the dependence on ||\1/m||oo can be replaced by ||^m||oo- 


Proof. We will restrict to the proof of the first statement, since the others are 
proved similarly. 

We will proceed in a Stampacchia iteration. Let p be a cut-off function with 
p = 1 on B{x, r/16) and p = 0 outside B(x,r/ 8) such that 0 < p < 1 and 
| V? 7 1 < 16. Let u := \x\ and for k > 0 set := max{r/?i — k, 0}. In addition set 
A{k) := suppufc. Then clearly A{k) C B(x,r/ 8), such that on A{k) the estimates 
from theorem 6.7 hold. 

The L 2 -bound for |y| from lemma 6.2 implies that 


k 2 \A{k)\< f p 2 u 2 dp, < f p 2 u 2 dp < C(r 1 , a, ||^m||oo, ||-^ E ||oo) • 

J A(k) JB(x,r/8) 

Therefore there exists ko = /co(|E|, ||\I/m||oo, \\K S ||oo, Co) < oo, such that \A(k)\ < 
Co for all k > ko. Here we want Cg to be the constant from lemma 6.8, to be 
able to apply the estimate from there for all functions with support in A(k), with 
k > ko. 

To proceed, let q > 2. Multiply the Simons identity, in the form (6.2) from the 
proof of proposition 6.5, by u q k and integrate. This yields 

[ —u g k uAu + uld'Vxl' 2 ~ |Vm| 2 ) dp 
J A(k) 

< c f u q k u A + \K\u q k u 3 + | L Rm s | u q k u 2 + u q k x * V(Q + x * S) dp . 

J A(k) 

Here c is a purely numerical constant. Partially integrate the Laplacian on the 
right hand side and the last term on the left hand side. This yields 


/ q7juu q k 1 \\7u\ 2 + qu 2 u q 1 (Vp, Vu) + u q k \ V%| 2 dp 

I A(k) 

<cf u g k u 4 + \K\u q k u 3 + | L Rm s | u q k u 2 + (u k \Vx I + u^u] Vu k \) (\Q\ +u|S'|)d/x. 

JA(k) 
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Note that the term f qrjuu 9 k ~ 1 \ Vu| 2 d/i on the left hand side controls f u q k \Vu\ 2 dfi. 
But before we use this estimate, we absorb the gradient terms into the first term 
on the left hand side using \Vu k \ < r/Vu + C(r~ 1 )u. Consider for example the 
term containing |Vx| 2; 

of ul\V X \(\Q\+u\S\)dfi< f ul\V X \ 2 dv + c[ u q k \Q\ 2 + u q k u 2 \S\ 2 d/i. 

J A(k) J A(k) J A(k) 

The other terms which contains Vn can be treated similarly, such that the 
resulting terms can be absorbed on the left. This yields an estimate of the form 


I A(k) 


u q \V u \ 2 dii<C(q,\\K*\U\\Ql 


-Rnr 


Ml 


I A(k) 


<7 — 2 

L k 


< 7-2 2 

a 


U 


d/i. 


( 6 . 8 ) 


Note that we used that u k < u and u < u 2 + 1 here to get rid of the extra terms. 
We begin estimating the terms on the right hand side of (6.8) using lemma 6.8. 
Rewrite and estimate the first term as follows: 


ul 2 u 6 d/i = 


(u k u 6 / q 


-2\q-2 


d/i 


I A(k) 


' A{k) 


< \A(k)\(c s q _ 2 / | V(u k u 6 / q ~ 2 )\ 2 + \Hu k u 6 / q ~ 2 \ 2 dfi 

V JA(k) 


q- 2/2 


(6.9) 


To estimate the first term on the right hand side compute on A(k), using u k /u < 1, 


|V(ii fc ii 6/,? - 2 )| < ii 6/,? - 2 |Vii fc | + ^ 2 ii 6/9 - 2 |Vii|^ < c(g,r- 1 )(u 6 / <? - 2 |Vu|+u^ +1 ). 

Observe, that if q is large enough, namely such that 2 + < p Q anc l J2_ + 2 < 

2 +po, then theorem 6.7 yields that 



V{u k u e/q ~ 2 ) 2 d/i < 


C(q ) [ u 12/q ~ 2 \Wu\ 2 + u^ +2 dfi < C(q). 
J A(k) 


Here, and for the remainder of the proof, C(q) denotes a constant that depends 
on q and, in addition to that, on all the quantities the constant in the statement 
of this theorem depends on. 

To address the second term in (6.9), recall that since 0 = 6 + = H + P, we have 
I|tf Hoc = HPHoo < 211^11^. Therefore 



< C(q) , 


where the last estimate also follows from theorem 6.7 if q is large enough. Sum¬ 
marizing these steps, we have 



u 6 d/i < C(q)\A(k)\. 
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A similar procedure for the remaining terms in ( 6 . 8 ) finally yields the estimate 



q k \Wu\ 2 < C(q)\A(k)\ , 


( 6 . 10 ) 


provided q > qo is large enough. Fix such a q > q 0 and let / = u k q ^ 2 . Since 
\Vf\ 2 <C{q,r- l ){ul\yu\ 2 + ulu 2 ), 
equation ( 6 . 10 ) and the above estimates imply that 



\Vf\ 2 dqi<C(q)\A(k)\. 


The Hoffman-Spruck-Sobolev inequality from lemma 6 . 8 , combined with theorem 
6.7, furthermore yields 


' A(k) 


f d/x = / uffi 2 djn < C(q)\A(k)\ 


' A(k) 


/ | V u | 2 + Hu 2 d/i 

'A(k) 


«3'd~2 

2 


<C(q)\A(k)\. 


Thus one further application of lemma 6.8 yields 


/ u q k +2 dn= [ f 2 dn < C(q)\A(k)\ 2 . 

J A(k) J A(k) 

Consider h > k > ko, then on A(h) we have that > h — k and therefore we 
derive the following iteration inequality 


\h — k\ q+2 \A(h)\ < f u q+2 dfi< [ u\ +2 dn < C{q)\A(k)\ 2 . 

J A(h) J A{k) 

The lemma of Stampacchia [Sta 66 , Lemma 4.1] now implies that \A{ko + d)\ = 0 
for 


d q+2 <C(q)\A(k 0 )\ <C(q)a 

In view of the definition of A(k) = supp ma x{tju — k , 0}, this yields in particular 
|x|(x)| < C(q). As x was arbitrary, the claim follows. □ 

Corollary 6.11. Let Sc M be an LM-stable MOTS with (r, a)-locally uniformly 
finite area. Then S satisfies the following estimates. 

sup |x| < C(a,r -1 , ||AT||oo, ||VA'||oo, || M Rm||oo, (inj(M, g))- 1 ) 
s 

In addition, for all x G S we have an L 2 -gradient estimate of the form 
I |Vx| 2 d/i < C(a,r~\ \\K\U \\VK\U I^RmlU). 

J B(x,r/8) 
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Proof. The statement to prove is that the constants only depend on the stated 
quantities. This is due to the following reasons. 

First, we use the above estimates using the representation of constants containing 
||Tm||oo an || M Rm|| 0O . As we have seen in remark 1, we can estimate 

|Tm| < c(|A'| 2 + |VA'| + | M Rm|), 

where c is a numerical constant. Second, since for all A", Y, Z e A(E) 

Q(X, Y, Z ) = L Rm(A, 1", n, Z) + L Rm(A, Y, is, Z), 
we can use the Gauss and Codazzi equations of the embedding M L to estimate 
\Q\ + | L Rm s | < c(|Af + | VAT| + | M Rm|). 

Third, obviously 

|AT S | 2 + \S\ 2 < \K\ 2 . 

Thus we see that all quantities are controlled by ||A'H^, ||VA'|| 0O and H^RmH^, 
where the oo-norms are computed on E. Note that the dependency on inj(M) 
comes from the fact that the constants Cq and cf in the Hoffman-Spruck-inequality 
only depend on || M Rm|| 00 and inj(M) -1 . □ 

We conclude with an estimate for the principal eigenfunction to Lm or A_. 

Theorem 6.12. Let E be an LM-stable MOTS. Let A > 0 be the principal eigen¬ 
value of Lm and f > 0 its corresponding eigenfunction. They satisfy the estimates 

A|E| + i f /- 2 |V/| 2 d / n<4vr+ f |S| 2 d,u- f * M dn 

«/ S «/ s s 

and 

I |V 2 /| 2 d^<C(|E|,||if|U,||Vii'|U|| ,, Rm|| 00 ,mj(M, 9 r 1 )£/ 2 + |V/| 2 d^ 

+ A 2 / / 2 d/j. 

The same estimates hold for L_-stable MOTS when f and A are the principal 
eigenfunction and eigenvalue of A_ instead, then has to be replaced by in 
the first estimate. 

Proof. The first estimate follows from a computation similar to the proof of lemma 
5.6. 

The second estimate then follows from the first by using the identity 

/ | V 2 /| 2 d/i = f (A /) 2 + s Rc(V/, V/) d/x. 

As As 
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To estimate the terms on the right hand side, note that 

-A/ = A/ - 2 S(Vf) - /(div 5 - Hxf - |S| 2 + f Sc - * M ) 
and as E is two-dimensional 

s Rc(V/,V/) = fSc|V/| 2 . 

In view of the Gauss equation for E C M and the bounds for x, we find the 
claimed estimate. □ 

Corollary 6.13. If E is an LM-stable MOTS, then the principal eigenfunction 
f > 0 to Lm which is normalized such that ||/||<x> = 1 satisfies the estimate 

J / 2 +|V/| 2 +|V 2 /| 2 d M < C(|E|, |E|-‘, || A"||oo, l| VA'Ilco, irRmlU.mjfM.j)- 1 ) 

The same estimate holds for L_-stable MOTS, when f is the principal eigenfunc¬ 
tion to L_ instead, and the constant depends on |E| ; |E| _1 ; \\K Hoc, I^RmHoo, and 
inj {M,g)- 1 ). 

Proof. Since ||/||oo = 1, we have f E f 2 dp < |E|. Then since f~ 2 > 1, the hrst 
estimate from the previous theorem implies 

J |V/| 2 < C(|E|, ||A|U, || VA'IU, irRmlU). 

Since A 2 / v / 2 < A 2 1E| < G|E| _1 the above estimates combined with the previous 
theorem imply the claim. □ 

Remark 6.14. A local version of these estimates can also be derived from local 
area bounds, like the curvature estimates before. In the subsequent application, 
however, we will not use this more general form. 


7 Local area bounds 

This section is devoted to derive the area bounds needed for the curvature esti¬ 
mates in the previous section. 

The following theorem is analogous to Pogorelov’s estimate for stable minimal 
surfaces [Pog81]. We will modify the proof of Colding and Minicozzi given in 
[CM02], 

Theorem 7.1. Let ( M,g ) be a Riemannian 3-manifold with curvature bounded 
above || M Rm|| 00 < C. 

Let E C M be an immersed surface with bounded mean curvature ||i7||oo A C. 
Assume that there exist a > 0 and a constant Z > 0 such that for all functions 
T) G C£°(E) we have 

— f s Sc ? 7 2 dp, < f (2 — a) | V ?7 1 2 + Zg 2 dp. (7.1) 

Jt, J t, 
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Then there exists r 0 = r 0 (a 1 ,Z, ||//’|| 0 o, U^RmUoo) such that for all r < r 0 , the 
area |E fl B^(x,r)\ is bounded by 

|E fl Bx(x,r) | < ^-r 2 . 


Proof. Fix an arbitrary x E T,. By the Gauss equation for E we know that 


s Sc = M Sc - 2 m Rc{u, u) + H 2 - |H| 2 < C'dl^Hoo, | M Rm|| 00 ). 


Hence by the Rauch comparison theorems (cf. [CE75, Section 1.10]), there is 
a radius 0 < rq = ?'i(||-H'||oo> || M Rm||oo) such that E has no conjugate points in 
B%(x, ri). Hence the pull-back 7 of the metric of E to the disc D ri := B( 0, ri) in 
T x E is regular and satishes (7.1). 

Denote D s = R(0,s) the disk of radius 0 < s < r\ in D ri and T s = dD s the 
boundary. Note that D s is a topological disk and is a single circle. Furthermore 
the area of D s with respect to 7 is bigger than \B%(x, s) |, for s < r\. 

In the stability inequality, we set = r)(s) = max{l ——, 0}, where 0 < r < r\ 
will be chosen below. Denote by K(s) = f D Seal dju. Then K\s) = f r Seal dl, 
where dl is the line element of r s induced by 7. Hence, by the co-area formula 
and partial integration 

K'(s)rj 2 (s) ds — I K(s)(r] 2 (s)y ds. 

J 0 


— / Seal if d/a = — 


1 o T 


Let l(s) = length(r. s ). By the formula for the variation of arc length and the 
Gauss-Bonnet formula we End l'(s) =27r — K(s), as is a circle, on which the 
geodesic curvature integrates to 2ti. Thus we compute, using (rj 2 (s)y = — 2 (1 —-), 
(?/ 2 (s))" = ^2 and the co-area formula 


Seal if dfi = 2n / (ri 2 (s)y ds + / /(s)(?/ 2 (s))" ds 


1 D r 


Jo 
2i r 

r 


O 1 _ “) ds + fZ 


l(s) ds = —27 t H—-| D r 


Furthermore, compute 


IV 7 I 2 d/x— — l(s) ds = —\D r \, 


1 Dr 


and estimate 


Zi] 2 d/i < Z\D r 


1 D r 


Hence equation (7.1) implies that 


-27t + \ \D r \ < \D r \ + Z\D r 


Thus 


a|D,.| < 27rr 2 + Zr 2 \D r \. 

Choose r 2 = min{^,r 2 } and absorb the error term. This yields the claim. □ 



The above theorem yields the local area bound for stable MOTS needed for the 
curvature estimates. 


Corollary 7.2. Let E be an LM-stable MOTS. Then there exists r 0 > 0 depending 
only on \\K ||oo, || V/l ||oo? and || M Rc|| 00 , such that for every i6E and r < r 0 

|E 0 B^(x, r)| < 67rr 2 . 

Proof. To see that theorem 7.1 is applicable on an L M -stable MOTS, recall equa¬ 
tion (5.1). As |y| 2 > \\A\ 2 — 4|A' S | 2 and by the Gauss equation 

|A| 2 = M Sc - s Sc - 2 m Rc(z/, u) + H 2 , 

we find from (5.1), taking the scalar curvature term to the left, that 

— T / S S c ?7 2 d/i< / 2\ V?71 2 + Zif dp 

2 J s Jt, 

where Z = Z^K^, HVATIU, I^RclU). As on a MOTS U^IU = ||P|U 
21| 77 s || oo, theorem 7.1 gives the desired bounds. 


Theorem 6.10 and corollary 6.11 imply the following estimates. 

Corollary 7.3. Let E be an LM-stable MOTS, then 

IMIco < C(Mco, IIVATIU J| M Rm|| co ,(inj(M, 9 ))- 1 ). 

Furthermore, there exists 0 < r — f(||A"||co. IIVA’m, || w Rm||co) such thxit for all 
x G E 


' B(x,r) 


IVx| 2 < c (Halloo, IIVATIU, || M Rm|| 00 ). 


8 Applications 


The main application of the curvature estimates proved in this paper is the fol¬ 
lowing compactness property of stable MOTS. 

Theorem 8.1. Let (g n ,K n ) be a sequence of initial data sets on a manifold M. 
Let (g, K) be another initial dataset on M such that 

|| M Rm|| 00 < C , 

Halloo + irv^iu < c, 

inj (M,g) > Cr 1 , 

for some constant C. Assume that 

g n -> g in Cf oc (M,g) and, 

K n —>• K m ClJM,g). 
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Furthermore, let £ n C M be a sequence of immersed surfaces which are stable 
marginally outer trapped with respect to (g n , K n ) and have an accumulation point 
in M. In addition, assume that the £ n have uniformly locally finite area, that is, 
for all x G M there exists 0 < r — r(x) and a = a(x) < oo such that 

|£ n fl BM tn (x,r )| < ar 2 uniformly in n, (8.1) 

where B Mtn (x,r) denotes the ball in M around x with radius r. 

Then a subsequence of the £ n converges to a smooth immersed surface £ locally 
in the sense of C 1,a graphs. £ is a MOTS with respect to (g, K). //£ is compact, 
then it is also stable. □ 

Proof. By the estimates in corollary 6.11, the above assumptions, even without 
(8.1), are sufficient to imply that the shears \ n , and thus the second fundamental 
forms A n , of the £ n , with respect to the metric g n are uniformly bounded 

14*1 < C. 

As g n is eventually C' 2 -close to g , this bound translates to a bound for the second 
fundamental forms A n of £ n with respect to the metric g. 

In the sequel Bm(x, s ) denotes an extrinsic ball in (M, g). As the geometry of the 
(M, g n ) is uniformly bounded, the uniform curvature bound implies in particular, 
that there exists a radius s such that for every x G £ n , the connected component 
of £ n fl B m (x, s ) containing x can be written as graph of a function u x n over T x ll n , 
where the function u * is uniformly bounded in C 2 . Without loss of generality, we 
can assume that s < r, where r is from equation (8.1). 

Now let x € M be arbitrary. From the previous fact we conclude that each 
connected component of £ n fl Bm(x, s ), which intersects £ n fl Bm(x, s/ 2) contains 
a uniform amount of area. In view of the local area bound (8.1), we conclude 
that there are only finitely many such components. Furthermore, the maximal 
number of those components is uniform in n. 

Hence, in each ball Bm(x,s/ 2), we can extract a convergent subsequence of the 
£ n , such that £ n fl B M (x,s ) converges in C l,a to a smooth surface EA As the 
M tn can be covered by countably many such balls, a diagonal argument yields a 
convergent subsequence of the £„ and a limit surface E, which is immersed (cf. 
remark 8.2). Note that since the £ n have an accumulation point in ( M,g ), the 
limit E is non-empty. 

Furthermore, C 1,a convergence yields that E is C 1,a and satisfies a weak version 
of the equation 6 + = 0. In view of standard regularity theory for prescribed mean 
curvature equations, we find that E is in fact smooth (cf. [GT98]). 

If E is compact, we can cover E with finitely many balls B(x,s/ 8). As before 
we know that locally the £ n converge to E in C 1,a . Since we also have local 
W 1 /2 -bounds on y we can furthermore assume that £„ —> E in W 2,p for a fixed, 
large p, which will be selected below. 
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From this we can conclude that, the metrics of the E n converge to the metric of 
E in C a fl W 1,p . We can pull back the metrics of the E n to E and call them y n . 
The metric on E will be denoted by 7. Then define the operators L n as the pull 
backs of the operator L M on E n to E. Let f n be the principal eigenfunctions of L n 
with eigenvalues A n and normalize such that ||/n||oo = 1- Since the area of the E n 
is eventually bounded below by half of the area of E, theorem 6.12 implies that 
0 < A n < C , where C = C(C, \\K ||oo, ||V/i ||oo, || M Hm|| 00 ). Thus we can assume 
that the X n converge to some A with 0 < A < C. 

By corollary 6.13, the VF 2,2 -norm of the f n taken with respect to the metrics y n 
is uniformly bounded. Recall that the difference of the Hessian of / with respect 
to 7" and 7 is of the form 


(V 


2 

7n 


7 )/ = 


■(r 7 „-r 7 )*tf/ 


where T 7 and T 7n denote the connection coefficients of 7 and j n . Furthermore 
V/ is bounded in any L p and by W 1,p convergence of the metrics T 7n — T 7 —> 0 in 
L p . Thus we fold that also ||/„||yc 2 . 2 < (7, where the norm is taken with respect 
to the metric 7 on E. Hence we can assume that f n —> f in W 1,p . The Sobolev 
embedding W 1,p < —> C°, implies that f > 0, and ||/||oo = 1, so / ^ 0. 

The next step is to take the equation L n f n = A n f n to the limit. Since f n —> / 
only in W 1,p , we have to use the weak version of this equation, namely that for 
all <j) G C^E) 


/ 7 n ( df n )id(j>j + B l n {df n )i(l) + C n f(j) d/i = X n / f n <j> d fi , 

J s J s 

where B n and C n are the coefficients of the operator L n . By the W 2,p -convergence 
of the surfaces, we find that y n converges to 7 in W 1,p , and B l n and C n converge 
in L p to the coefficients B l and C of Lm on E. Thus, since f n converges in W 1,p 
to /, we can choose p large enough to infer that the limit of the above integrals 
converges to the corresponding integral on E, that is / satisfies 


/ (V/, V0) + (B, Vf)cp + Cf<j)dp = X f(j) dp . 

.7s J s 

Thus / is a weak eigenfunction of Lm on E. Elliptic regularity implies that / is 
smooth and satisfies L M f = A/. Since A > 0 and / > 0, / ^ 0, we conclude that 
E is stable. □ 

Remark 8.2. If the limit surface is not compact, then it still follows that it is 
“symmetrized” stable in the sense that inequality (5.1) holds for all test functions 
rj with compact support. 

Remark 8.3. (i) If the surfaces E n are embedded, one would assume that the 
limit E is embedded as well. However this is not necessarily the case. This is 
due to the fact that at a point p where E touches itself, the equation 9 + = 0 is 
satisfied with respect to the outward normal. At p the normals corresponding to 
these two sheets point into opposite directions at p. If we flip one of the normals, 
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to make them point into the same direction to apply the maximum principle in 
a graphical situation, the equation for one sheet will remain 9 + = 0 , but for the 
other it will change to 6~ =0. Hence, one cannot compare the two sheets. 
However, if E also satisfies 6~ < 0, then the maximum principle implies that the 
set S of touching points is open. By continuity S is closed, and S ^ E, as the E n 
are embedded and have bounded curvature. Hence 5 = 0 and E is embedded. 

(ii) If the assumption of uniformly locally finite area does not hold, but the 
surfaces E n are embedded, the limit still exists in the sense of laminations. Here, 
the limit is a lamination, for which the leaves are “symmetrized” stable MOTS 
(cf. remark 8.2). Convergence is in the sense of laminations in the class C a , 
with convergence of the leaves in C 1,a , for any 0 < a < 1. For a proof of this 
statement, we refer to [CM04, Appendix B], in the reference this is stated for 
minimal laminations, but the modification to the MOTS case is straightforward. 
For the sake of completeness we state the definition of a lamination. A lamination 
C C M is a closed set, which is a disjoint union of complete connected smooth 
surfaces, called leaves. Furthermore there are coordinate charts for M, ip : V C 
M —> R 3 , V a neighborhood of some point x G M, such that the image of each 
leaf L of C is contained in a set of the form R 2 x t, where t G I, and / is a closed 
subset of R. A sequence of laminations C n is said to converge to a lamination 
C if the coordinate charts converge and C is the set of accumulation points of 
the C n . 

(iii) There are examples of compact three dimensional manifolds which contain 
sequences of compact stable minimal surfaces of fixed genus and unbounded area 
[CM99], [Dea03]. Thus assumption (8.1) does not follow immediately from stan¬ 
dard theory, even when the surfaces E n are confined to a compact region. 
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